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Colloidal particles absorbed at the interface of a liquid droplet arrange into unique packings during 
slow evaporation (Manoharan et al. Science 301 483-487). We present a numerical and theoretical 
analysis of the packing selection problem. The selection of a unique packing arises almost entirely 
from geometrical constraints during the drying. 



Manoharan et al. 1] recently presented an ingenuous 
method for fabricating clusters of small particles into pre- 
cise configurations. Polystyrene spheres (diameter 844 
nm) were dispersed in a toluene- water emulsion with each 
oil droplet containing a low number N of spheres. The 
toluene was then preferentially evaporated, forcing the 
particles to come together into compact clusters. Sur- 
prisingly, the final particle packings were unique: the ob- 
served packings ioi N < 11 closely correspond to those 
previously identified 2] as minimizing the second mo- 
ment of the particle distribution, M = ||r^ — ro|p, 
where tq is the center of mass of the cluster. 

The fact that such a simple process leads to preci- 
sion assembly at the submicron scale points to exciting 
possibilities for controlling the assembly of more gen- 
eral objects 3]. The goal of this Letter is to under- 
stand the physical principles underlying the observations 
of Manoharan et al. Why are the final packings 

unique? (Why) do they minimize the second moment? 
What physical parameters do these results depend upon? 
We first present numerical simulations of hard spheres on 
an evaporating liquid droplet for a wide range of liquid- 
solid contact angles: for each contact angle the simula- 
tions reproduce the final packings of ^ . We then demon- 
strate that the uniqueness of the packings, as well as their 
connection to minimal moment structures, can be under- 
stood from purely geometrical considerations. The argu- 
ments suggest a methodology for creating new packings, 
which we confirm through numerical simulations. 

Numerical Simulations. The experiments suggest the 
following theoretical problem: for a given liquid volume, 
the particle configuration is determined from minimizing 




(a) 



(b) 



(c) 



(d) 



FIG. 1: Schematic two-dimensional representation of the dry- 
ing experiment: (a) particle configuration when the droplet 
volume is above the critical volume; (b) critical packing; (c) 
rearrangement below the critical packing with both capillary 
and contact forces acting on each particle; (d) final packing. 



the total surface energy 

U^=1D I dS^-iDP I d^ + 7p / d5, (1) 

JD JdP Jp 

while respecting excluded volume constraints between 
the particles. Here DP and P refer to the droplet 
surface, droplet-particle interface and particle surface re- 
spectively (see Figure nj. 

Numerical simulations of this problem are performed 
using Brakke's Surface Evolver 4], a program which de- 
termines the equilibrium configuration of deformable sur- 
faces given the definition of an energy. The colloidal 
spheres are modelled as liquid droplets with high sur- 
face tension, typically 1 to 2 orders of magnitude larger 
than the main droplet, in order to penalize non-spherical 
deformations of their shape. Interfacial tension between 
the droplet and the particles are chosen appropriately in 
order to satisfy Young's law at the solid-liquid contact 
line, 7p = jdcosO -\- Jdp^ where 6 is the equilibrium 
contact angle. Non-interpenetrability is enforced with 
an excluded volume repulsion energy Ur acting between 
the centers of the spheres 12j; Ur dominates when at 
least two spheres overlap by one percent. 

The particles are initially positioned randomly on the 
droplet. The droplet volume is then slightly decreased 
(by one percent or less) and the particles rearrange to 
a new equilibrium, minimum of U = -\- Ur. This 
procedure replaces the evaporation dynamics by a series 
of equilibrium problems and therefore mimics the low 
evaporation rate limit of the experiments. 




FIG. 2: Comparison of the experimentally observed packings 
(left) with those obtained by numerical simulations (right) as 
a function of the number of spherical particles. 
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FIG. 3: Numerical evolution of the second moment as a func- 
tion of the droplet volume for = 9 (both in units of the par- 
ticle radius). The final second moments are 25.899 (theory) 
and 25.946 (simulations). Insert: packings observed during 
the drying process when M = 25.946, 27.016, 35.014, 45.058. 
Solid line: minimal second moment; dashed line: second mo- 
ment at the critical packing {M = 27). 



We find that the packings obtained numerically are 
unique and agree with those obtained by Manoharan et 
at. m, over the range of contact angles tested (10° < < 
170°) and initial conditions. A comparison of the final 
computational and experimental packings is illustrated in 
Figure 121 In all cases, the final second moment obtained 
numerically differs by less than 0.5 % from the minimum 
moment packings of 2]. The agreement verifies that 
the experimental packing problem corresponds to the en- 
ergy minimization posed above and that other particle- 
particle interactions (electrostatic, van der Waals...) are 
not essential to generate the final packings. 

There is one noteworthy difference between the simula- 
tions and the experiments. In the simulations, the cluster 
of spheres evolves smoothly and continuously throughout 
the drying process (c.f. Fig. 01 for 9 particles). In con- 
trast, experiments [l| showed a discontinuous transition 
to the final packing, with an abrupt change in the second 
moment at a critical volume. The experimental disconti- 
nuities must arise from experimental features not present 
in the simulations, probably contact angle hysteresis. 

The Critical Volume. We now study theoretically the 
packing selection problem. For sufficiently large liquid 
volumes, the minimum energy solution is a spherical 
droplet with noninteracting force- free particles (Fig. ^). 
There exists a critical volume Vc below which the droplet 
cannot remain spherical (Fig. ^p); below this critical vol- 
ume each particle is acted upon by capillary forces. At 
the critical volume Vc^ there is a critical packing of par- 
ticles, which can be characterized as follows. 

The interactions between particles on the surface of a 
sphere are equivalent to the steric interactions between 
the "cone of influence" of each particle, defined as the 
cone originating from the droplet center and tangent to 
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TABLE I: Characteristics of critical packings of spherical col- 
loidal particles as a function of their number N: number of 
contacts (ric), number of independent forces (n/) and number 
of admissible modes of rearrangement (rim)- 

the particle (Fig. ^p). These are also equivalent to the 
interactions between the intersection of these cones with 
the droplet, which are circles. Consequently, packing 
spherical particles at the critical volume is equivalent to 
packing circles on a sphere, a mathematical problem with 
rich history H, @, 0, and for which numerical solutions 
have been proposed up to = O(IOO). Usually, the cir- 
cle packings for a given N are unique, with two types 
of exceptions. For N = 5, 19, 20, 23, 26, 28, 29 . . . the cir- 
cle packing has continuous degrees of freedom (where at 
least one circle is free to "rattle") 8]; for A^ = 15 two dif- 
ferent configurations (15a and 15b) lead to the densest 
packings with equal surface density 8]. 

Hence, provided that there are no kinetic traps, identi- 
cal particles will arrange themselves into the circle pack- 
ing at the critical volume. Fig. jSlconfirms this conclusion 
in numerical simulations for A" = 9. Note that changing 
the contact angle of a particle changes the size of its 
circle of influence and the droplet radius by the same 
amount, so the circle packing problem is unchanged and 
the critical packing of particles at the critical volume is 
independent of wetting characteristics; note however that 
the critical droplet volume is contact angle dependent. 

How do the particles rearrange when V < Vc? The en- 
ergy minimization problem suggests that we must find 
the particle configuration which minimizes Uy: under 
inter-penetrability and contact angle constraints. Ow- 
ing to the complexity of solving for the liquid surface of 
constant mean curvature, this at first appears to be an ex- 
traordinarily difficult theoretical problem. However, we 
have found that in fact the constraints associated with 
packing of particles are sufficient to uniquely determine 
the initial rearrangements of the particles. 

Let us suppose the droplet volume is reduced by a small 
amount SV <C Vc. Deviations of the droplet interface 
from spherical lead to capillary forces (F^) on each par- 
ticle. Since every particle must be in force equilibrium, 
these forces must be balanced by contact forces (fj^) be- 
tween the particles (Figure QJ): 

F» + E = 0' (2) 
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3 4 4 4 

4 6 6 6 

5 9.333 9.333 9.333 

6 12 12 12 

7 16.683 17.100 16.683 

8 21.157 21.657 21.157 



9 25.899 25.899 25.899 

10 31.828 31.828 31.828 

11 37.835 37.929 37.835 

12 43.416 43.416 42.816 

13 51.316 52.690 47.701 

14 59.225 60.279 54.878 



TABLE II: Final second moment of the model (Mm) com- 
pared with the final second moment in the experiment 1] if 
the particles are assumed to be perfectly spherical (A^exp) 
and with the minimum second moment (A^2)- 



where C{i) denotes the set of particles in contact with 
sphere i. 

Let us characterize the number of ways the parti- 
cles can rearrange to accommodate this change in vol- 
ume. Each particle has three degrees of freedom and the 
droplet has one (the value of its mean curvature or pres- 
sure), so there are 3A^ + 1 degrees of freedom. The con- 
straints are of three types: (a) solid body rotation does 
not modify the packing (3 constraints); (b) the particles 
cannot overlap, {uc constraints, where Uc is the number 
of contacts at the critical packing); (c) forces have to 
balance (Eqn. [^J. 

Equation Q implies nontrivial constraints. Indeed, 
suppose we try to solve (0) for the ric contact forces 



(e. 



Gj, where e,; is the unit vector 



directed from the droplet center to the center of the par- 
ticle). Equation Q has SA/" components, and ric ^ 2N 
unknowns (Table I). Consequently, solutions exist only 
if compatibility relations are satisfied between the F^. 
Since capillary forces depend on the position of the par- 
ticles, these equilibrium considerations constrain the re- 
arrangement of the particles. Geometrical constraints 
have also played an important role in understanding force 
propagation in granular packings 9]. 

Close to the critical volume, the capillary forces are 
given by = F^ei. From the N scalar forces F^, equa- 
tion ^ shows that only a subset n/ can be chosen inde- 
pendently, i.e. equilibrium of each particle lead to N — rif 
additional constraints; n/ is found by computing the rank 
of the compatibility matrix in (O 13] . 

The total number of admissible modes of rearrange- 
ment rim for the colloidal particles is found by subtract- 
ing the number of constraints from the number of degrees 
of freedom: we find rim = 2N — 2-\-nf — nc. For a given 
rim is entirely determined by the geometry of the critical 
packing. The results are displayed in Table ^ When- 
ever the circle packing is unique, we find that rim = 1- 
When rim > 1, we find that the number of modes is al- 
ways correlated with the presence of rattlers. If a total 
number of continuous degrees of freedom rid exist in the 
circle packing {nd=2 or 1 per rattler depending on if it 
completely free, as in = 19 or constrained in a slot, 
as in A/" = 5), we always find that rim = 1 + ^d- Since 
rid indicates that there exists rid force-free surface 




FIG. 4: Comparison for N = 9 between the packing given by 
(a) experiments, (b) simulations and (c) model. 



modes for the packing, we obtain therefore that there 
exists a unique mode of rearrangement for all particles 
which are non-rattlers. In experiments, the degeneracy 
in the circle packing problem is chosen by additional in- 
formation: for example in Q, because the particles are 
charged only on the side exposed to the water, there is a 
weak dipolar repulsive force between the particles which 
breaks the degeneracy [l^ . 

This result implies that there is only a single set of 
{Fi} that is consistent with all the constraints. This 
mode is independent of surface energies, and depends 
only on geometry. Capillarity does enter into the problem 
in relating the force Fi to the displacement Sri of the i*^ 
sphere. If radius of curvature of the droplet changes from 
R to R -\- SR then volume conservation of the droplet 
implies that Svi and 5R are related through A Svi + 
{4:71 R^ — NA)SR = 0, where A is the wetted area of the 
particles. The capillary force Fi is then given by 



_ ^ / r Acos 3 



Sri 



(3) 



where a is the dry angle of the particle on the critical 
packing, the equilibrium contact angle, and /3 = a — 0. 
This formula is asymptotically valid in the limit 6V 
so that deviations from a spherical cap droplet are small. 

The above results apply just below Vc. However, the 
general principle can be iteratively applied below the crit- 
ical volume: starting from the spherical packing, we de- 
crease the droplet volume by small increments and, as- 
suming equation ^ continues to hold, we calculate the 
corresponding incremental particle rearrangement con- 
sistent with all the constraints. At each step in the it- 
eration it is necessary to recompute the mode {Fi}; for 
every packing there is a unique choice that is consistent 
with the constraints. This process iterates until the final 
equilibrium configuration is reached. During the process 
the packing changes substantially from the initial disk 
packing- typically multiple new contacts are added. 

The results of the model are displayed in Table UTI and 
in Figure 21 The model reproduces accurately the final 
experimental packings for N < 6 and 9 < A/" < 14; in 
particular, the non-convex packing for A^ = 11 is well 
predicted by the model. There are small differences in 
the final packings for A^ = 7 and 8, likely arising from 
deviations of the capillary force-particle displacement re- 
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(a) Cb) (c) 

FIG. 5: Final configuration of 9 spfieres witfi different wetting 
conditions: on the droplet, ni spheres have a contact angle 
of 160° and n2 = 9 — ni spheres have a contact angle of 20°. 
(a) cases {711,712} = {2,7} and {4,5}, M = 27.706; (b) case 
{711,712} = {6,3}, M = 29.780; (c) case {711,712} = {8,1}, 
M = 30.754 (units of the particle radius). 

lationship from the linear law Overall, the model 
leads to a dimensionless error in second moment with 
the experimental packings of 1.6% versus 2.6% for the 
minimum second moment criterion. It should be empha- 
sized that the computational cost of the new algorithm is 
orders of magnitude slower than that of a full simulation. 
To our surprise the final packing can be computed quite 
accurately without ever knowing the shape of the liquid 
surface during the packing process! 

Finally, we remark on the minimal moment criterion 
itself: our results suggest that the drying influences the 
final packings only through (a) enforcing the initial disk 
packing at the critical volume; and (b) through equa- 
tion ((2|), relating the Fi to the particle displacements. 
It is probably not coincidental that (up to prefactors) 
equation Q is similar to the force-displacement relation 
Fi = —diAi ^ c\5ri + C2 ^^i? minimal moment 



criterion. We believe that this explains the similarity be- 
tween the observed structures and those minimizing the 
second moment. 

The calculations presented herein suggest that the 
unique packings observed by Manoharan et. al. arise 
because (i) the initial circular packing is unique for the 
regime they explored, except for A/' = 15; and (ii) the sub- 
sequent evolution of the particles is so highly constrained 
that there is only one final packing that is consistent with 
the constraints. 

This suggests that the only way to generate different 
packings is to modify the circle packing at the critical 
volume. This can be easily modified by choosing par- 
ticles with differing sizes or wettabilities (thus creating 
circles of different sizes). We have run simulations where 
particles on a given droplet possess different contact an- 
gles. Figure shows three cases of A/" = 9 particles with 
contact angles of either 160°or 20°. The three packings 
and their second moment differ significantly from Fig. ^ 

In summary, we have presented a numerical and theo- 
retical study of the packing selection problem of Manoha- 
ran et al. 1] . The selection of a unique packing was found 
to arise almost entirely from geometrical constraints dur- 
ing the drying process and a procedure was proposed to 
generate different packings. This could be extended to 
characterize all admissible evaporation-driven packings 
of colloidal spheres. 
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